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Abstract 

Many Lagrangians of physical theories can be expressed as eigen- 
values of certain, relatively simple, matrices involving Dirac gamma 
matrices. We give concrete examples for Lagrangian corresponding 
to a point particle coupled to electromagnetic field, electrodynamics, 
nonabelian gauge theories, extended objects and gravity. We also dis- 
cuss (in case of a point particle) what are the implications of the least 
action principle applied to matrix Lagrangians. 

1 Motivation 

A nice property of Dirac matrices in 4 dimensions is that a large class 
of Lagrangians appearing in classical physics can be understood as 
eigenvalues of some matrices involving them. As for the Lagrangian 
for a point particle this observation is almost trivial considering the 
context in which Dirac introduced gamma matrices - it also works 
in arbitrary spacetime dimension. In the following sections we would 
like to note some less trivial observations, in particular that the La- 
grangians of 4D gauge theories and general relativity as well as ex- 
tended objects can also be represented as eigenvalues of certain ma- 
trices. Lastly we consider the variation procedure for actions that are 
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matrices (not eigenvalues of matrices) and discuss conditions under 
which one recovers classical equations of motion. 

The fact that Dirac matrices appear so naturally at the level of 
classical theories (although they originate from Dirac's attempt to 
formulate relativistically invariant quantum theory) is somewhat in- 
triguing. If not purely accidental and formal it must be a consequence 
of some underlying physical reason. For the letter a good explanation 
would a claim that the spin (and hence the Dirac/spin matrices) is the 
property of space-time itself which is then inherited by particles. Un- 
der these circumstances it would be natural to expect that relativistic 
field theories involve (in a certain way) Dirac matrices already at the 
classical level. 

2 Lagrangians as eigenvalues 
2.1 Point particles 

For a point particle with mass m, charge q and coordinates x^, coupled 
to the electromagnetic field A^x), in curved spacetime given by the 
metric g^ u (x), one simply takes 

P := m^{x)x>* + lqA^x)^, {^(x), 7l/ (x)} = 2g^(x). (1) 

Here 7 ^(x) are x dependent Dirac matrices. One possible choice is 
to use the vielbein e"(x), g^, v {x) = e a Ax)e\,{x)r] a b where 7] a i, is the 
Minkowski metric (we use the many minus convention when needed) 
and to take 7 M (x) = e^(x)7 a where 7a are usual (x independent) Dirac 
matrices {<y a ,"fb} = ^Vab- 

Because the square of the mass term of P is just m 2 g fll/ (x)x fl x u 
and since the potential term commutes with the mass term we con- 
clude that there are two doubly degenerated eigenvalues of P and the 
spectrum is 

spec{P) = P± = ±myj g^WxKx" + qA^x)x». 

Hence the Lagrangian for a point-particle can be obtained by taking 
the smallest eigenvalue, P_, of P so that the least action principle 
gives the correct equations of motion. 
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2.2 Gauge theories 

A similar result may be obtained for a Lagrangian of electrodynamics 
in flat space in 4 dimensions. One takes 

E := l -{! ab F ah f = -\F ab F ab l - \i^F ah *F ab (2) 

where r y ab := py a , 7 6 ]/2, F ab := d\ a Au is the field strength and the dual 
tensor is *F ab := ^e abcd F c d where e abcd is completely antisymmetric 
with e 0123 = 1. The second equality in ([2} is unique for 4-dimensional 
spacetimes and relies on the following identity for gamma matrices 

{ 7 afe , 7 cd } = -2(7] ac rj bd - V ad rj bc )l - 2i^e abcd . (3) 

The eigenvalues of E are again double degenerated and its spectrum 
is 

spec(E) = E± = -^F ab F ab ± \iF ah *F ab . (4) 

(where we use the fact that there are two doubly degenerated eigen- 
values of a matrix a + 7 5 6: a± b). The second term in E± is a total 
derivative and hence does not affect the equations of motions there- 
fore to obtain the formulation of electrodynamics one can take E+ or 
eM Before going further let us note that by equation (j2J) it is not 
possible to couple a pseudo-scalar field (the axion) to the F*F term, 
leaving the FF term unchanged. Therefore from this perspective the 
coupling to the axion field does not appear here quite naturally. 
For non-abelian gauge theories we may take 

ill 

E Y M • = -|(7 G afe) = ~^ G ab G 1 - G ab G 

where A is the adjoint index of the gauge group and G^ b is the A'th 
component of the field strength (reducing to F ab if the group is U(l)). 
Here the conclusions are similar to with one difference that the 
G*G term must be kept due to the presence of instantons. Therefore 
the 9 angle is fixed and equal to 9 = 8tt 2 (using the conventions where 
the 9 term is given by-J^G^ *G ab A ). 

One can also take * F a b instead of F a b in the definition of E obtaining equivalent result. 
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2.3 General Relativity 

For the Lagrangian of General Relativity we find that a good choice 
is 

G := \R lu , pa 'f v {xyf a {x) = R1- 2i 7 5 *R = R1. (5) 

where 7 A " y (x) = e^e^ ah and *R is a dual to the Riemann tensor. To 
prove that ([5]) holds one again uses ([3]). Moreover, in 4 dimensions 
that dual *R is zero due to the Bianchi identities 

*R ■- e^ pa R^ pa = 8(^0123 + i?023i + #0312) = 

therefore matrix G is actually a unit matrix proportional to the Ricci 
scalar. 

2.4 Extended objects 

For extended objects (with p spatial dimensions) the task is more 
tricky since the corresponding Lagrangians are more complicated (com- 
pared to a point-like particle) . Their Lagrangians are given by [U El [3] 

L = -AyJ(-l)PdetG aP , G a p = d a X^X v g^(X) 

and the action is obtained by integrating L over the world-volume 
parametrized by a a , a = 0, . . . ,p (for strings, p = 1, for membranes, 
p = 2). Here X^ are the coordinates of the object and g flll {X) is the 
metric of spacetime. Dimensional constant A is the tension. Using 
observations made in [3] we find that a good choice of the matrix 
Lagrangian is 

P s trin 9 = ^v(X){X»,X' / }, (6) 

{X»,X U } = e afi d a X^dfiX v 

for a string and 

^membrane = ~^'y[ivp{X){X P ' , X ,X P }, (7) 

{X P ,X U ,X P } = e a ^d a X^d l3 X u d 1 X p 

for a membrane, where ^ V P{X) = e^e^ abc , j abc = \{^ v y + cycl). 
We note a natural appearance of the Poisson and Nambu brackets in 
((6|) and ([7]) respectively. For higher dimensional objects p > 2 (which 
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would require the embedding spacetime to have more dimensions then 
4) one simply takes j^ T yJ tn ... f , p+1 {Xt M , . . .,X^}. 

Let us now argue why the eigenvalues of matrices © and ([7]) give 
the corresponding Lagrangians. Using identity ([3]) one immediately 
verifies that the square of P s tHng is equal to — detG Q/ g, a,/3 = 0,1. 
To prove the same result for a membrane it is useful to use another 
identity 



ry^ 7 a'6' c '| _ 2e abcd e a ' b ' c ' d 'r lddl 


1 = -2 


^aa' 


rf> 


Yjbc' 






if*' 


cb' 




which follows form the fact that in 


4 dimensions 


rylbc 


= ie° 



Now it is easy to see that the square of P mem brane is equal to det G a p , 
a, j3 = 0,1,2. 



3 Action principle for matrices? 

Let us apply the least action principle directly to a matrix Lagrangian 
(not its eigenvalue) of a neutral particle in curved space 



S = m J 7 ^(x)x^dr. (8) 

In the case of a flat space the matrix Lagrangian in S is a total deriva- 
tive hence 5S = for all trajectories - which is of course inconsistent 
with the geodesic equation. 

In curved space the conclusion is different. Varying S with respect 
to and integrating by parts we obtain 



5S= I {x»d ul ^{x)-%)5x u dT (9) 

hence the equations of motion are x^dy^^x) = % which can be also 
written as 

= 0, x=dtfy v \{x). (10) 

However the only solution of this equation is x^ = which can be 
seen by expanding the matrix fn V in terms of 7 a 's as 

ffiu = f^ula, f^v '■= d[v e v \ 
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therefore equation (|10p implies x^f® = for every o hence x^ = 0. 
This shows that the direct application of the least action principle to 
matrices via (JSj) is inconsistent with the geodesic equation. 

In the remaining part of this section we argue that the geodesic 
equation can nevertheless be obtained for certain geometries if the 
equations of motion (jlUp are modified as follows 

x^Uv = C„. (11) 

Here C„ is some constant, nonzero matrix and can be obtained by 
evaluating the l.h.s. of (jlip at some time tq. Therefore matrix C u is 
completely determined by initial conditions of the trajectory x^(t) via 
C u = x^ f fMU \ T=TQ . This modification could be argued from equation © 
if we generalize the variation procedure by saying that the integrand in 
is at most a traceless (but nor necessarily 0) matrix. The matrix C v 
can be expanded in terms of 7 a 's and therefore satisfies this condition. 

Let us now find solutions of (jlip . Ideally, we would like to find the 
inverse of f^ v so that equation (fTT|) could be inverted. This can be 
done by noting that in 4 dimensions there exist a unique identity for 
antisymmetric rank-2 tensors and their duals, namely [5] 

/£, *f up a = ~(fua *r a K, a fixed, 

where *f^ u " is a dual tensor *f^ u a := e^ CT /^/2. This identity 
implies that 

iu», *r p } = -2^(/7), iff) ■■= a *r b va b , (12) 

where *f» u is a dual matrix ■=* a la . Applying (p|) to (fTTI) 
we find that 

^ = -^y{^*/n (13) 

where we assumed that the geometry under consideration is such that 
(/*/) 7^ 0- That assumption implies that the metric g^vix) is non 
standard (i.e. not diagonal). To see this explicitly let us assume that 
the metric is given by g^ u {x) = r] tiv f^{x) (n fixed) hence e° = 5^f^(x). 
In this case it follows that (/*/) is 

(/*/) = [-(5o-i/ / ,)(9 2 -3/,)+(9o-2^)(5i- 3 /,)-(5o-3/ M )(ai- 2 /,)]r/^ 

(where 9 a _b := d a — d^) which is 0. Therefore equation (fl~3|) applies to 
very special geometries. 
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Equation (|13p is an autonomic one i.e. it gives the formula for 
the 4- velocities x^ in terms of x^'s. To obtain the solutions for 
one would have to integrate (|13p once, which however can be done 
explicitly only when the metric is specified. 

Let us now show that x^ satisfying (|lip satisfies also the geodesic 
equation. To see this it is useful to differentiate (modified) equations 
of motion (jlip w.r.t. r. We find that 

^ ffiu "b X^X dcrf^v 0. 

Applying (TL21) to the above equation we obtain 

x" + Y^x»x° = 0, % := *r a d { Jl )v Va b . 

Therefore we arrive at the geodesic equation provided T^a coincide 
with Christoffel symbols r^ CT . 
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